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Abstract 

Short-term variability of X-ray spectra has been reported for several Active Galac- 
tic Nuclei. Significant X-ray flux variations are observed within time scales down to 
10 4 — 10 5 seconds. We discuss short variability time scales in the framework of the 
X-ray flare model, which assumes the release of a large hard X-ray flux above a small 
portion of the accretion disk. The resulting observed X-ray spectrum is composed 
of the primary radiation as well as a reprocessed/Compton reflection component 
that we model performing numerical radiative transfer simulations. We conduct 
Monte-Carlo simulations of large flare distributions across the disk including rel- 
ativistic corrections of the observed radiation and discuss general implications on 
AGN spectra and their variability. 
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1 Introduction 



The X-ray emission of Active Galactic Nuclei (AGN) can be explained by two 
components: a primary radiation emerging from a hot, optically thin medium 
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and a reprocessed component coming from colder, optically thick medium. The 
existence of these two media can be explaine d by th e flare model, which was 
first suggested by Galeev. Rosner. fc Vaiana ( 1979f ) and developed in many 
subsequent papers for both galactic black holes and AGN. The AGN flares 
are assumed to be driven similarly to solar flares: magnetic reconnection cre- 
ates a compact source above the accretion disk, the source emits the primary 
spectrum partly reaching the observer directly and partly illuminating the 
disk atmosphere, the atmosphere then emits the reprocessed spectrum. 



In this paper we present the first modeling results for the reprocessed radi- 
ation coming from the irradiated area, the so-called hot-spot, underneath a 
compact flare source as indicated in Fig. 1. We take into account details of the 
flare geometry including the dependence on the incident angle of the primary 
radiation. 
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Fig. 1. Flare geometry 

We further conduct Monte-Carlo simulations of random flare distributions 
across the accretion disk in order to construct rms X-ray spectra. These simu- 
lations also contain relativistic corrections due to the curved space-time in the 
vicinity of the black hole. We show general trends for several parameters such 
as the radial distribution of hot-spots and of the incident radiation flux over 
the disk. 



2 Modeling 



2. 1 The reprocessed spectra across the hot-spot 



We perfor m detailed radiative transfer us i ng th e codes TITAN and NOAR de - 
scribed in Dumont. Abrassart. fc Collin ( 20001 ) and in Dumont et al. (2003). 
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TITAN conducts radiative transfer inside optically thick media, determin- 
ing the temperature and ionization structure. NOAR is a Monte-Carlo code 
treating the Compton processes. We consider a plane-parallel section of the 
accretion disk atmosphere at 18 R g (R g = M being the mass of the 
black hole) from the disk center. The hydrostatic equi librium of such an at- 
mosphere is calculated following iRozariska et al. (1999) a ssuming a 10 8 solar 



mass black hole and an Eddington accretion rate of 0.001. Collin et al. (2003) 



have shown that a change of the hydrostatic equilibrium due to the onset of 
the flare happens slowly compared to the adjustment of the temperature and 
ionization profile inside the medium. We therefore assume that the density 
structure of the medium does not change for a significant fraction of the total 
flare duration. 



For the incident flux F inc we presume a power law with F inc oc E~ a , a = 0.9, 
between 0.1 eV and 100 keV. The ratio F inc / Fdisk of the incident flux to the 
flux coming from the disk is fixed at 144, hence the flare by far dominates over 
the disk emission. 



It can be seen from Fig. 1 that the incident angle and the ionization parameter 
£ = IzlFW f the primary radiation (with Fi nc being the incident flux and tin 



being the hydrogen density at the disk surface) depend on the position within 
the hot-spot. This aspect has not yet been considered in p revious flare models 
where the disk illu mination was defined as semi-isotro pic dCzernv et al. L 2004 ; 
Collin et al. . 20031 ) or with a constant incident angle ( Navakshin fc Kallman . 



200ll) . 



For reasons of energy equipartition the distance h of the flare source above 
the disk surface should be simil ar to the pressure scale h eight of the disk 
( 
( 



Galeev. Rosner. fc Vaiana . 1979j ). Applying equation (2.8) of lShakura fc Sunvaev 
19731 ) to our case we derive h ~ 1.9 x 10 cm ~ 0.013 R g . We define a half- 



opening angle Go = 60° for the flare cone, which, for an isotropically emitting 
source, corresponds to 50% of the incident flux reaching the disk being repro- 
cessed. Then, the hot-spot radius of h x tan 60° ~ 0.023 R g verifies our implicit 
assumption that the hydrostatic equilibrium computed at 18 R g is sufficiently 
accurate across the whole hot-spot. 



The observed spectrum is constructed as the integrated spectra of successive 
concentric rings at incident angles with 0° < < ©o- The spectra can be 
evaluated at various disk viewing angles i (measured from the disk symmetry 
axis). Details of the calculation method will be published in a following article 
(Goosmann et al., in preparation). 
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2.2 The rms-variability of random flare distributions 



To investigate the consequences of this more detailed modeling of the hot-spot 
we perform Monte-Carlo simulations of many flares d i stribu ted over the ac- 
cretion disk. The procedure is outlined in lCzernv et al.l (2004J ). We assume the 
hot-spot emission to be independent of the flare position. Radial probability 
distributions p(Ri, <pi) for the appearance of a flare can be defined as 



p(Ri, fa) = p{Ri)p{<f>i 



p7rad+2 TDlrad+2 



1 

27T" 



Herein R4 and <j>i are polar coordinates, Ri n and R ou t denote the inner and the 
outer disk radius. For 7 TO <z = the hot-spot distribution is uniform while for 
positive (negative) values of 7 the number of flares rises in the outer (inner) 
disk part. We also enable a radial dependence of the hot-spot incident flux: 

(2) 



Since the dissipation within a Keplerian disk scales with r~ 3 we expect j3 ra d 
of the order of 3. Finally, the radial dependence of the flare life-time is chosen 
to scale with the Keplerian time scale. 



The average number of hot-spots existing on the disk is denoted by n mean} 
the mean X-ray luminosity of the primary source by Lx = 8 x 10 43 erg s _1 . 
The rms-spectra are computed from n seq = 50 consecutive model observations, 
each of an integration time T in t,,„ = 6156 s. Relativistic effects are i ncluded by 



sec 



Dovciak. Karas. fc Yaqoobl 2004 . for a code 



application of the code KY 
description). For the following computations we assume a Kerr black hole of 
10 7 solar masses with spin parameter a = 0.95. The inner disk radius is at 
1.8 R, 

to be i = 30°. 



g , the outer one at 50 R g . The viewing angle of a distant observer is taken 



3 Results 



3. 1 The reprocessed spectra across the hot-spot 



In this paper we only examine the pure reprocessed component (e.g. without 
dilution by the primary radiation) for various positions within the hot-spot, 
e.g. for various incident angles of the irradiating flux. In Fig. 2 we show the 



4 



reprocessed flux EFe versus photon energy E at several incident angles G. 
The results are shown for a face-on view, as assumed in Seyfert-1 galaxies, 
and for an intermediate viewing angle, as supposed for Seyfert-2 galaxies. For 
both viewing angles spectral variations are found in the soft X-rays across 
the hot-spot. These changes are more significant at higher viewing angles and 
affect the whole spectrum. Note that the dependence on the incident angle of 
the primary spectrum is opposite for the two viewing angles shown. At face- 
on view the flux decreases with rising incident angle, while at intermediate 
viewing angles it rises. 




log (E/[l keV] ) log( E/[l keV] ) 

Fig. 2. Reprocessed spectra for different incident angles of the irradiating hard 
X-rays. Dotted line: cos 9 = 0.95 (9 = 18°). Solid line: cos 9 = 0.65 (9 = 49°). 
Dashed line: cos 9 = 0.35 (9 = 70°). Left: face-on view. Right: intermediate viewing 
angle at 62° from the normal. The inlets show a zoom between 0.4 and 0.8 keV. 

The spectra emitted by the whole hot-spot are shown in Fig. 3. For increasing 
viewing angles the spectra decrease systematically, as does the soft X-ray slope 
and the ratio of the Compton hump to the soft X-ray spectrum. 
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Fig. 3. Left: integrated spectra of a hot-spot with half-opening angle of 60°. The 
viewing angles increases from top (i = 13°) to bottom (i = 89°). Right: flux ratios 
versus viewing angle for such a hot-spot. Dashed line: F[5.5keV]/F[4keV] (soft X-ray 
slope). Solid line: F[28keV]/F[4keV] (hard/soft X-ray slope). 
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3.2 The rms-variability of random flare distributions 

In Fig. 4 (left) we show the dependence of the rms spectra on the radial 
distribution of the incident radiation flux. Again, we do not consider any 
dilution by the primary radiation of the flare. We assume that the hot-spots 
at outer radii are dimmer, but if the decrease of the luminosity is fast, the rms 
spectrum has a much more complicated structure since it depends strongly on 
a few innermost brightest hot-spots. Computing various random realizations of 
the process without any variations in the model parameters, we see statistical 
differences between them. If the hot-spot's surface brightness decreases more 
slowly the role of the outer hot-spots increases, causing a major drop in the 
overall variability level. 




E [keV] E [kcV] 

Fig. 4. Left: Dependence of the rms/mean spectrum on f3 ra d (other model parame- 
ters: j ra d = —1.5, n mean = 300). Right: Dependence of the rms/mean spectrum on 
7 rarf (other model parameters: (3 rad = 3, n mean = 100). 

In Fig. 4 (right) the rms spectra for several radial hot-spot distributions are 
shown. The value of j ra d = corresponds to a random but uniform coverage 
of the disk surface by the hot-spots. If the flares and hot-spots appear more 
frequently in the outer part of the disk (7 ra d > 0), the overall variability 
decreases but a variable red wing is still present. If the flares appear relatively 
more frequently in the inner part of the disk (7 ra d) < 0, a clear dip is seen at 
the position of the unshifted iron line. 



4 Discussion and Conclusions 

Modeling the reprocessed radiation of the hot-spot reveals that the pure re- 
flection component varies with the position in the hot-spot. This variation is 
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stronger for intermediate (edge-on) viewing angles than for a face-on viewing 
direction. The effect is caused by the £ parameter varying with the distance 
from the hot-spot center. For outer parts of the hot-spot £ is lower, it will 
thus cause a colder medium and a thinner hot surface of the atmosphere. The 
temperature of the medium determines its absorption properties, while the 
hot surface will influence, by electron scattering, the angular dependence of 
the re-emitted flux. The incident angle of the radiation also has an impact on 
the spectral shape as the Compton scattering phase function is not isotropic. 

The integrated spectra show a strong dependence on the viewing direction. 
This effect is partly explained by the size of the projected surface area of the 
hot-spot seen at different inclinations. However, there also is a dependence of 
the overall spectral slope on the viewing direction that could be used to put 
observational constraints on the disk inclination. 

The modeling of rms spectra for distributions of many flares across the disk 
reveals that flares at the inner part of the disk dominate the variability. Either 
a strong concentration of flares at the inner disk part or a relatively higher 
luminosity of inner disk flares increases the level of variability. Due to the light 
bending, an orbiting flare in the central region of the disk is effectively seen at 
various local viewing angles, which increases spectral variations. Although the 
overall level of variability is only weakly dependent on the energy, there are 
variations of the rms spectra for the energy range of the Ka line. This effect 
is due to the strong gravity close to the Kerr black hole. It could give a hint 
to t he lack of variability in t he iron Ka line observed in MCG -6-30-15 (see 
e.g. Fabian fe Vaughan . 20031 and references therein). 



A basic assumption of our modeling is that the density structure of the medium 
stays constant after the onset of the flare. Hence, we give results for flares 
having a short duration compared to the dynamical time scale. On longer 
time scales the medium exp ands to obtain a new hydrostatic equilibrium (see 



Czernv fe Goosmanrl 12004 for details) and the spectral properties change. 
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